THE DERIVE - NEWSLETTER #63 

ISSN 1990-7079 


THE BULLETIN OF THE 



1 Letter of the Editor 

2 Editorial - Preview 

3 User Forum 

6 Halloween in St. Polten 
Tania Koller 

7 Thick Lines and Colour Pots in DERIVE 
Bernhard Kutzler 

15 TI-Nspire - The next Generation? 

Josef Bohm 

25 Beyond Polynomial Regression (DERIVE and CAS-TI) 

Peter Schofield & Josef Bohm 
32 Contour Plots (& NSOLUTIONS improved) 

Bill Wilburn 

40 Special CAS-Guest: MATHEMATICA 
A Non Polynomial System of Equations 


October 2006 









D-N-L#63 


Bookshelf & Information 


D-N-L#63 


[1] Einfuhrung in die Matrix-Algebra / Mit Anwendungen in der Statistik, Karsten Schmidt & 
Gotz Trenkler, Springer 

[2] Wiskunde voor hoger onderwijs door middel van Derive 1 & 2, Peter van der Velden et al.. 
Academic Service - SDU, The Hague, NL 

[3] Experimente und Visualisierung, mit Beitragen von G. Schmidt, E. Lehrmann, J. Bohm . , 

MU Der Mathematikunterricht, Jg 52, Heft 4-2006, Erhard Friedrich Verlag, Seelze 

[4] Dynamisch Geometrie entdecken, Elektronische Arbeitsblatter fur die 7.-10. Klasse, 

H.-J. Elschenbroich u. Gunter Seebach , coTec-verlag 

[4] and many other valuable books (German and English) are available from 

http://shop.bk-teachware.com 


Recommended websites: 

For your information: 

http://www-fourier.uif-qrenoble.fr/~parisse/enqlish.html 

Daniel Legendre 



Nils Hahnfeld was busy as ever ... 

GREAT NEWS: Calculus Made Easy version 8.0 is now available. 

Calculus Made Easy v8.0 has become the ultimate educational tool to study Calculus. 

No other software provides detailed Step by Step Solutions to Integration (using U-Substitution 
(choose u yourself or let CME do it for you!), Integration by Parts, by Partial Fractions, Power Rule, 
Expand/Rewrite, Arctan(x) and Arcsin(x) Integrals, etc) and Differentiation (Chain-, Power- , Product- 
and Quotient Rule) problems. 

Just enter the given function and learn how the correct rule is applied - step by step - until the final 
answer is derived. 

No other software covers Calculus I (Calc AB) , II (Calc BC) , III (multivariable), Differential Equations, 
Vector Calculus and Differential Geometry!! 

To view some examples, check here: www.ti89.com/documentation.htm 
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letter or the evitor 



Dear DUG-Members, 

with some delay and with many I present 
DNL#63. I must admit that I intended to skip 
DNL#63 and intended to have a double issue 
63/64 in December but then I prepared a paper 
for a seminar - the nonlinear regression - and 
from one moment to the other I wanted to put 
it into a smaller than usual DNL#63. As you 
can see, it did not luck. This issue is again a 
full publication. Many mails from and to Peter 
Schofield, Karsten Schmidt, Bernhard Kutzler, 
Tania Koller and Bill Wilburn were exchanged 
and I am happy that the 44 pages are again full 
of materials. 

In this DNL you might have your first meeting 
with the new TI-CAS device and software TI- 
Nspire. Bernhard Kutzler gives a short intro- 
duction. If anybody of you has some experi- 
ence with Nspire, he or she are friendly wel- 
come to share it with us. The developers an- 
nounce that as many DERIVE features as pos- 
sible will be implemented into Nspire. I hope 
that they will have enough inspiration to do it! 

Several DUG-members asked if we could not 
open our publication for other computer alge- 
bra systems. Of course, we can and we start 
with our first special CAS guest, 
MATHEMATICA. Bill Wilburn sent a respec- 
tive contribution. It might be nice for students 
to reproduce calculation and graph with 
DERIVE (or the TI!). So if you have any arti- 
cle in Maple, MuPAD, Mathematica, Auto- 
graph (which is not a CAS, but a very fine 
math program) etc, then please send it to me. 

I complained for years that DERIVE is unable 
to plot in different line thicknesses. Tania 
Koller shows us to overcome this deficiency. 
The “trick” is so easy and I am very disap- 
pointed that it was me who found it. Tania's 
students are very eager using thick lines and 


also combining DERIVE and a graphic pro- 
gram - which is also very novel. (I am quite 
sure that many of you had never tried this). 
Tania sent a collection of Halloween-Ghosts 
which were created by applying polynomial 
functions. They are all descendants of Paul 
Drijver's ghost - Hi Paul, do you remember 
Krems 1993, where you presented your ghost? 
Your ghost feels very well and has survived in 
our regions (maybe because of the excellent 
wines of Krems and surroundings?) 

Admire the creativity of Tania's students on 
pages 6 and 14. 

One contribution is redemption of my promise 
of DNL#59 to produce contour plots in 
DERIVE. This is in a tight connection with 
implicit plots and so we can have implicit plots 
with thick lines. Thanks Peter Schofield for 
your great cooperation. It was just yesterday 
that Peter sent his improved tool for numeri- 
cally solving equations which is used for cal- 
culating the level curves. (More about this in 
DNL#64.) 

Some interesting requests are remaining for the 
next DNL (Slider bars and Arrows, a very fine 
tool to collect coordinates of 2D-points in the 
Algebra window, a TI-bug?, . . .). 

I received a fine article from G P Speck, New 
Zealand, who rejoined the DUG, welcome 
again. 

Finally I'd like to invite you to order the Pro- 
ceedings CD of the DES-TIME Conference. It 
is full of materials (papers, programs, files, 
pictures, ...) 

Best regards until December 



Download all DNL-DERIVE- and Tl-flles from 

http : / / www . austromath . at /dug/ 
http : / / www . derive-europe . com/ support . asp?dug 
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I 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with contents of 44 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE 
and the 77-89 /92/ Titanium! Voyage 200 as 
well as to create a group to discuss the 
possibilities of new methodical and didac- 
tical manners in teaching mathematics. 

As many of the DERIVE Users are also 
using the CAS-TIs the DNL tries to com- 
bine the applications of these modem tech- 
nologies. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: December 2006 

Deadline 3 1 November 2006 


Preview: Contributions waiting to be published 

Two Stage Least Squares, M. R. Phillips, USA 
Some simulations of Random Experiments, J. Bohm, AUT 
Wonderful World of Pedal Curves, J. Bohm 
Another Task for End Examination, J. Lechner, AUT 
Tools for 3D-Problems, P. Luke-Rosendahl, GER 
ANOVA with DERIVE & Tl, M. R. Phillips, USA 
Financial Mathematics 4, M. R. Phillips 
Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon & Co, J. Bohm 

Challenges from Fermat, Bj. Felsager, DEN 
Are all Bodies falling equally fast, J. Lechner 
Modelling Traffic Density, Th. Himmelbauer, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Mathematics and Design, H. Weller, GER 

Diophantine Polynomials, D. E. McDougall, Canada 

Challenger Matrix Problems, G P Speck, New Zealand 

Precise Recurring Decimal Notation, P. Schofield, UK 

and Setif, FRA; Vermeylen, BEL; Leinbach, USA; Koller, AUT; Baumann, GER; 
Keunecke, GER and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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Three mails from Bill Wilburn: 

I am glad to hear that the DNL will continue on. A while back, when I heard that Dr. Rich and Theresa 
were no longer with Derive, I suspected that this might happen. For a long while, I ? ve been interested 
in Mathematica. So, I downloaded the trial version. I liked it. My wife got me Mathematica 5.2.0 Pro- 
fessional for Father's Day. I'm still a novice with it. It has so many features. It's amazing, and you can 
do so much with it. I personally like all of the CAS systems. I've had no experience with Maple, but it 
appears to be a fine CAS system as well. I hope that you can get the latest version of Mathematica. I 
know that a person of your caliber certainly deserves it. 

Sincerely, 

Bill Wilburn 

PS. I'm looking forward to more issues of the DNL. 


Hi Josef: 

Here's an example of a Mathematica Notebook on solving 2 sets of parametric equations. 

I used DERIVE 6.10 as a word file to show the Mathematica Notebook's contents. 

I suppose that you could put it into a DNL if you wish. If so, please feel free to edit it as you like. 

Have a nice day, 

Bill Wilburn 


Hi Josef: 

I ran across this site and thought that you might be interested: 

http : / /curvebank .calstatela.edu/ famous curves/ famous . htm 

Regards, 

Bill Wilburn 

You can find Bill's example on page 38 together with the respective DERIVE treatment , Josef 


Manuel Garcia Mateos 


garcia_mateos@WEB . DE 


Hello Derivians, 

I tried to show that the sequence (l+l/n) A n <3 for all natural numbers by defining n as a positive inte- 
ger and solving the inequality. 

Derive does not resolve the inequality for n. Derive produces the following output: 


#1 : n :e Integer (0 , «) 

#2: SOLVE 


((* * t)" dtr) * LNW 


Is there a way to show that the sequence is bounded and monotonuously increasing without using Ber- 
noulli-inequality or other transformations? 


Regard the sequence e(n):=3*(-0.7) A n (n positive Integer) The sequence is bounded and converges to 
zero, but: solve(e(n)<5,n) leads to false! 

Is there a way to show that alternating sequences are bounded? 


Greetings, 

Manuel 


p 


4 


Derive - and CAS-TI-User Forum 


D-N-L#63 


G P Speck, Wanganui, New Zealand 

Dear Josef: 

1 1 was good to receive your email and pleasing to hear that you remember Wanganui fondly. 

For various reasons I’m having difficulty sending technical material via email - so for this one time I 
am sending a disc with the Derive V.5 file eh_0(Challenger).dfw containing my paper Challenger 
Matrix Problems which you viewed previously. With this disc you can, of course, modify the format 
however you wish for sending to DUG members. Not being familiar with Derive V.6 until I receive 
my copy, I hope that this will be no problem. 

In the meantime I’ll try to get my difficulties on margins, matrices, email formats, and other printing 
matters sorted out. 

Keep up the good work with Derive - we’ll hope that most of the good features of Derive will be 
maintained in NSpire. 

Sincerely, 

G P Speck 

Many thanks to G P. The diskette arrived in Wurmla and its contents will make a fine article 
fora new ACDC-column. Josef 

Peter Schofield, UK 

Hello Josef, 

Thank you for your news letter informing me of the end of TPs interest in Derive. I'm quite upset 
about it, but I suppose it was coming from the time they decided not to employ Albert anymore. I have 
strong reservations about "Nspire" and what it is intending to achieve. I fear that many of the wonder- 
ful aspects of Derive will be lost. 

Anyway, I am attaching a short article - together with corresponding Derive and TI-89 files - on how a 
precise recurring decimal notation can be set up in Derive 6 and on the TI-89. 1 hope this might be 
suitable for the DUG? 

All the best 
Peter 

Peter's contribution is on my stack and will be published in the next DA//.. 

Peter Witthinrich, Germany 

Lieber Josef Bohm, 

schade, dass die Entwicklung von Derive eingestellt wird. Derive ist in meinem Umfeld eher akzep- 
tiert worden als die CAS-Rechner, ich habe es im Unterricht deutlich haufiger als die Kleinrechner 
eingesetzt. 

Hoffentlich iibersieht TI nicht, dass so gut wie alle Schulen ohnehin einen Computerraum haben, sich 
aber nicht alle Schulen zusatzlich noch weitere Rechner anschaffen konnen. 

AuBerdem gibt es inzwischen eine Anzahl weiterer Computeralgebrasysteme, auch wenn Derive aus 
meiner Sicht besonders (oder besser als die anderen) fur den Schulbereich geeignet ist. Manche Lehr- 
krafte steigen vielleicht einfach um, zumal einige CAS kostenlos sind oder die Grundversion (die in 
vielen Fallen fur den Unterricht ausreicht) kostenfrei ist. 
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Meistens hat bisher der Elternverein die kleinen Rechner bezahlt. Jetzt werde ich den Eltern klarma- 
chen miissen, dass nach der noch nicht so lange zuriick liegenden Anschaffung der Voyage200 schon 
wieder ein neuer Typ besorgt werden soil. 

Zur DUG denke ich, dass sie doch ohnehin bereits auf CAS-Rechner erweitert war. So sehe ich keinen 
Grand, sie nicht bestehen zu lassen. Falls erwiinscht konnte sie spater immer noch erweitert oder um- 
benannt werden. 


Mit besten GriiBen 
Peter Witthinrich 


Lisa Townsley, USA, sent a short answer wrt my report on the "(non) Future" of DERIVE in my DUG- 
Information. I am very happy these three sentences. Many thanks on behalf of the DUG, Lisa. We will 
go on. 

Josef, 

This was a nice, sensitive note to send. I have shared your words with my department. The 
best part about using Derive is the excellent colleagues all over the globe. 

I sent Karsten Schmidt my paper about the "nonlinear" regression (page 25) and he sent back an ex- 
tended and interesting answer. I try to translate his mail (JB.dfw and the pdf-file can be downloaded). 

Dear Josef, 

in the meanwhile I dealt with your mail from, in particular with your first example ( n = 10). (page 25). 

The extremal value problem (P(x)=a*l+6*sin(x)+c*cos(x) is a linear regression with 3 exogene vari- 
ables (including a pseudo variable for the intercept). Even if sin(x) and cos(x) are created by a nonlin- 
ear transformation, the regression model remains linear. 

This can easily be observed by inspecting your matrix M, which corresponds to matrix X ? in chapter 
12 of our textbook (published by Springer, for an example see page 218) 1 . 1 reprodraced the procedure 
in DERIVE using the “usual” notation (file JB.dfw). #3 creates the regressor matrix X from x.; KQ- 
estimation with X and y leads exactly to your result. #7 and #8 allow a 2D-scatter plot together with 
the regression line. (Comment: I believe that KQ is "kleinste Quadrate” = "least squares”.) 

The 2. attachment is the pdf-SPSS output. 

1. page: Scatterplot and "linear” regression; 

2. page: nonlinear regression with model equation a+6*sin(x)+c*cos(x), i.e. the model from above 
with transformed variables, which leads to the same results; 

3. page: nonlinear regression with model equation a*sin(b*x+c)+d, which is equivalent to the SinReg- 
treatment of the TI-92+/Voyage 200, which gives the same result as in your paper.also entsprechend 
dem Ansatz SinReg des TI-92, es kommt auch das raus wie in deinem Text. 

Under "Residual”/” Sum of Squares” you can fin in SPSS the respective SSE - sum of squared errors - 
(twice 1.772 and once 1.645), which shows that the (true) nonlinear regression fits better. 

By the way, this is a fine demonstration of GauB-Markow-Theorem which says that the KQ-estimator 
(LS-estimator) is the best (i.e. minimizing variance) of all linear unbiased estimation function, because 
the better estimator is a NONlinear one! Best regards - Karsten 


see our book shelf. 


i 
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Hallo, it's Halloween in St. Polten! 
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Tania Roller s students like tftick lines, colours, 

inatiieiuaiics and Halloween 
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Thick Lines and Colour Pots in DERIVE 

Tania Koller, Austria 

My colleague Tania gave a great presentation in Dresden (on the Proceedings 
CD). Among many other “motivativing activities” for students she showed her “Cin- 
deralla's Shoe” modelled by cubic splines. I have seen this earlier, but look at her 
“DERIVE picture” of the shoe: thick lines and a bright pink, which is not among the 
DERIVE shading colours. I asked how she produced thick lines in DERIVE which 
has been a long requested feature of DERIVE, and how she brought the colour on 
the shoe. Here is her answer with some additional comments. Josef 




4 2 

x - 13 *x + 30 

#1: 

f(x) := 2 


30 

*2: 

TABLE(f (x) , x, -5, 5, 0.01) 

*3: 

TABLE(f ' (x) , x, -5, 5, 0.01) 

*4: 

f"00 


Set Option > Display > Points Connected and 
Large and Plot #2 (don't forget to set Options 
> Approximate before plotting. 

Then Set Option > Display > Points con- 
nected and Medium and Plot #3. 

Finally plot #4. Choose colours as appropri- 
ate. 



If you intend to use this nice tool more often then you might use a function with fixed parame- 
ters which can be changed if needed. 
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In function thick()variable and the increment of the table are preset with x and 0.01 . In most 
cases these defaults will be fine - and you need only entering function and boundaries: 

#5: thick(u, mi, ma, v_ := x, inc := 0.01) := TABLE(u, v_, mi, ma, inc) 

#6: thick(f(x), -3, 3) 



#7: thick(20.t + 1, -0.5, 0.5, t) 



I enter variable t and see that the steep- 
ness of the graph needs a smaller incre- 
ment. 


*S: thick(20*t + 1, -0.5, 0.5, t, 0.001) 



Now it looks fine. 

There is no Anti-Aliasing-algorithm - so set the 
points as tight as possible. 


The slider bar is one of my favourite features of DERIVE 6 (and as I know, it is also very of- 
ten used by Tania in classroom). So I am very happy that the "Thick- and Medium-Line-Trick" 
is wonderful supported by the slider bars. Sometimes it needs some moments to create the 
graph, but then the slider bars work without any delay. 


#9: 


thick(a*SIN(b*x) , - 4*jr, 4 *tt) 
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In parameter form we do not need this trick, because we can enter the size of points: 


*10: [2*SIN(3 *t) , 3 •C05(2*t)] 



Plotting in Size Medium we face 
problems because of the restric- 
tion to plot a maximum number of 
640 points. To obtain a smooth 
line one has to plot in intersections 
0 - 7t, jr - 27i, etc. 

#11: [2*SIN(5 *t) , 3 *C0S(2*t)] 

Or we try to overcome this restric- 
tion by defining a function applying 
the more or less unrestricting 
VECTOR-command: 



#12: thickpar(u, mi, ma, v_ := x, inc = 0.01) := VECTOR(u, v_, mi, ma, inc) 

#13: thickpar([2*SIN(5 *t) , 3.COS(2.t)], 0, 3 >n, t) 



Incited and inspired by Tania's method I wanted to do implicit plots in three line thicknesses, 
but this is not so easy to perform. One idea is to adapt the utility for plotting level curves (see 
another contribution in this DNL). I tried another way: Scanning the area in small x- and y- 
increments and plotting the points which fullfill the implicitly given equation of the curve: 
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impl(u, xst, xend, yst, yend, inc := 0.01, acc := 0.01, pts, i, j, val) := 
Prog 

pts := [] 

u := LHS(u) - RHS(u) 
i := xst 

Loop 

If i > xend 

#15: RETURN pts 

j := yst 
Loop 

If j > yend exit 

If ABS(SUBST(u, [x, y] , [i, j])) < acc 
pts = APPEI'D(pts, [[i, j]]) 
j :+ inc 
i :+ inc 


#17: 




X 

2.37 

y 

2.37 ” 

impl 



+ 


= 1, -4, 4, -4, 4, 0.03 


V. 

3.67 


2.88 

J 



(I'll come back to this “Superellipse” in the next DNL, Josef) 


I wanted to test my function and invented a strange implicit form and DERIVE gave immedi- 
ately its graph: 


3 - x*y 

x *e 

#18: = 1 

2 2 

x - y 


But we want the graph as a 
medium or as a thick line. 

Autograph is not a CAS but 
has a lot of great features, 
with a very powerful graphic. 

I plotted the expression with 
a thickness of 3 pts: 



And I took DPGraph to plot the left hand side surface together with the plane z = 1 for better 
illustration: 



# 19 : 



2 , 


- 1 . 


0 . 002 , 0.005 


v> X - y 



Another possibility is offered by adapting one or the other function of producing contour- or 
level plots (see the contour plot - contribution in this DNL). If you don’t mind calculation time 
then try the following and plot applying different point sizes: 


' 3 - 

x 


IMPPLOT 


x-y 


= 1 , - 1 , 2 , - 1 , 2 , 0 . 1 , 0 . 1 , 0.1 
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Now let's turn to Painting in DERIVE. Tania told me that one of her students (Eva Maria 
Korntheuer, 16 years) came across an Option in the Edit-Menu which was never mentioned 
in class before: Edit > Bitmap-Object. (Have you ever noticed this?). Just for fun she wanted 
to know what would happen. She converted the graph (eg a circle) into a Paintbrush Picture. 
And then with a double click on this picture she found herself in the Windows Tool Painbrush 
with all its features at her disposal. Many thanks Eva Maria and keep your interest and curi- 
osity. (She created the bear on page 14.) 



ij Derive 6 - [Algebra 1] 


E File 


D qS 


# 1 : 


Edit Insert Author Simplify 5o[ve Calculus Options Window Help 


Expression... 

Enter 

Text or Plot... 

Enter 

. ... 


Annotation . . . 

Links to OLE Objects 


Bitmap-Objekb 

X Cut 

Ctrl+X 

% Copy 

Ctrl+C 

8^ Paste 

Ctrl+V 

X Delete 

Del 

Undelete 

Ctrl+Z 

Select All 

Ctrl+A 

Mark and Copy... 

Ctrl+Shift+M 


m 


» s n. \ him 9 j- zn kj; 


Eearbeiten 

Offnen 


2 

y = 4 
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In the Graphics Program you can manipulate the graph - remaining in DERIVE - as you like 
and clicking in the DERIVE screen (eg on the expression) you see the coloured graph. Any 
double click on it leads back to Paintbrush for further treatment. 




Tania's students created a whole gallery of pictures using various functions. She sent a se- 
lection of pictures and I am very happy to present a part of this gallery in the DNL. 



I found the graph below in a mathematics book it is a chal- 
lenge for students to reproduce this pattern. 


r o.s-t 

spir := L - 0.5*e *SIN(t), 


0.5-t 

0. 5*e -COS(t). 


spir2 := 


rot(s) := 


t/2 t/2 

e • SIN(t) e -COS(t) 

2 ’ 2 

' COS(s) - SIN(s) ' 

. SIN(s) COS(s) . 


VECTOR 


^spir 


k • 7T 


k, 0, 17 , VECTOR 


spir2*rot^ j, k, 0, 1 7 J 


Plot with 0 < t < 7i . 


It is not so easy to find the mathematics behind. This is one way to produce a similar pattern 
by two families of curves. The graph was distorted in Paintbrush and then the tiles could be 
filled. 



You can change the DERIVE built in colour palette, but the shadings look pale because of 
the shading algorithm, the Paintbrush colours are really bright colours! And it would be very 
tiresome to shade the many facets of the pattern above ... and of the graphs on the next 
page! 
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Gallery of DERIVE-Graphs treated as Paintbrush-Pictures 

(produced by Tania's Students of the Handelsakademie St. Polten, Lower Austria) 



The dromedary (by artist Daniel Wailzer) is made up by a list of polynomials up to degree 10! 

You can find a lot of other students' products in Tania's great presentation given at 
DES-TIME 2006 in Dresden on the Conference CD. 


Unfortunately Tania did not unhide her secret how to bring 
colour on the screen of her Voyage 200!! Josef 


F3 Y FH 
^ace|Regraph 


ifriathlo 



[•EG HUTO 


FUNC 
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TI-Nspire - The next Generation? 

Bernhard Kutzler 

Linz, Austria 


At this year’s Derive conference Texas Instruments and Soft Warehouse Europe announced the dis- 
continuation of Derive. The reason is that TI will soon release their new CAS product TI-Nspire and 
in the future all development and marketing resources shall be focused on this new product. TI-Nspire 
carries the technology of the legendary TI-92 further into a new generation Math/CAS tool. 

This article is written for people familiar with Derive and presents some TI-Nspire features which are 
not available in Derive, hence can be regarded an enrichment by a Derive user. This text is based on a 
pre-release version of the TI-Nspire software (version 1.0.500). The same can be performed on a 
handheld, on which some of the screen images look slightly different. 

There are many features in Derive 6 which are not (yet) available in TI-Nspire. Most of these features 
are there, because the user community demanded them. These features resulted from the feedback of 
countless teachers who used earlier versions of Derive in their classrooms and reported their findings 
and suggestions to the Derive development team. Clearly, if you use such a feature in your classroom, 
you would miss it in TI-Nspire, in particular if you were among those who asked for it. 

There are exciting features in TI-Nspire which we had asked for in Derive, but our wishes were not 
(maybe could not be) granted. So the ideal would be a union of both products . . . 

In all areas of life we see technology leaps. From 1967 to 1972 NASA used Apollo rockets for their 
space flights with Apollo “versions” 1-17. Then they changed to the Space Shuttle program. Today we 
use cars with combustion engines which bum oil products. In a world with oil getting short, hence 
getting more and more expensive, we will soon need a “technology upgrade” to a next generation of 
motors using alternative fuels. Change is the essence of life! 

TI-Nspire is an exciting new product. Its modern design has a potential which goes far beyond what is 
practically possible in Derive. By offering, seamlessly integrated, computer algebra, dynamic geome- 
try, spreadsheet technology, word processing, data collection, and data exchange with a handheld, you 
now have, “under the same roof’, all technologies which are used today in mathematics teaching. 

My hope is that many Derive users will look at TI-Nspire, test it in the classroom, compare it to De- 
rive, and help with their feedback to shape TI-Nspire so that it will be not only Derive ’s declared suc- 
cessor, but will become an upgrade to Derive, hence a tme “next generation”. When comparing the 
two products we do need to bear in mind that Derive is an adult, while TI-Nspire is a new kid on the 
block, quickly growing up in today’s (software) world, certainly having a long and exciting future 
ahead. So please see the potential - and help turning as much as possible into reality. 
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TI-Nspire 


In the opening screen of TI-Nspire the preview pane on the left shows the overall structure of the 
document. A TI-Nspire document is a collection of “problems”. (An empty document starts with only 
one problem, named “Problem 1”.) Each problem is a collection of “pages”. Each page can be used as 
a whole or it can be divided into up to four sub pages. Each page (or sub page) can be connected to 
one of the following four applications: “Calculator”, “Graphs & Geometry”, “Lists & Spreadsheet”, 
and “Notes”. Right of the preview pane is the working area, which shows the selected “page”, hence 
offers an interactive environment for the respective application (or applications, if the page was di- 
vided into sub pages). 




The left screen image shows a page which was divided into three sub pages, but with no applications 
connected to them yet. The right screen image shows a second page added to the first page, with a 
menu open for connecting an application to it. 


Calculator / CAS 

The “Calculator” application provides the CAS and corresponds to the “Home” screen on a TI- 
89/92/92+ or Voyage200. 

When entering expressions one receives a lot of valuable support. Say, we want to input 
x + x + sin(x) . After typing x+x+si, the following is displayed on the screen: 

x+x+si 

Variables are displayed in italics and currently this would be interpreted as the sum of the three terms 
x, x, and si. As soon as we type the letter n, we see on the screen: 

x+x+sin 

TI-Nspire now recognizes the string sin as the name of a built-in function, hence does not display it 
in italics any longer. Adding the letter x leads to: 

x+x+si nx 


Now the third term, sinx, is all displayed in italics, because it is interpreted as the name of a four- 
letter variable. This is a clear “warning” to the user that this is not the intended input, which was the 
sine of x, simply because there need to be a pair of parentheses around the argument of the sine func- 
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tion 2 . We delete the last character, x, with the backspace key and add an opening parenthesis, (. 


x+x+sin 


x+x+sin( 


TI-Nspire automatically introduces a matching closing parenthesis with the cursor placed in between. 
Typing x completes the desired input: 


x+x+sin(x) 


The automatic insertion of a matching parenthesis is a great help to avoid “missing parenthesis” errors. 

Entering the input with the return key I produces a pair of expressions, comprising the entry and 
the answer, the latter being the (algebraically) simplified version of the entry. 


x+x+sin(x) 


sin(x)+2 


x 


Unlike Derive and like the TI-92/89 or Voyage200, TI-Nspire automatically simplifies every entry. 


For entering powers, derivatives, integrals, systems of equations, limits, sums, products, matrices, 
piecewise defined functions, and certain function logic constructs (such as if-then-else and loops) TI- 

Nspire offers two-dimensional templates. Say, we want to compute the derivative of x 2 ln(x) . We 
type x, then hit the keyboard’s power key 0 followed by the digit 2. 



The 2 is displayed as a raised-up exponent. The current position of the cursor indicates that the next 
character would be entered on the exponent level. The right arrow key exits the exponent level 

and brings the cursor down to the level of variable x. The tab key L“±lI would do the same. 


Since there is no ambiguity, we may simply type “In (x”, i.e. without a multiplication symbol (and 
also without the closing parenthesis). 



For computing the derivative we can choose the corresponding function (and template) from the Ex- 
pression Templates menu or the Calculus menu: 




*3[0[3 GGfflSBLK 

| Expression Templates'] 



2 TI-Nspire, like the TI-92/89/Voyage200, does not offer Derive’s Character Input Mode. 


pl8 


Bernhard Kutzler: TI-Nspire - The next Generation? 


D-N-L#63 


We open the Expression Templates menu, ... 


E EH® 10 S 







x+x+si 

x 2 -ln() 

l | 



. . . from which we choose the Derivative template. 



The red dots under and the vertical cursor left of the left-most (and lowest) dashed rectangle indicate 
that the next input goes into this position. Obviously this is for the variable. We type x followed by 

the right arrow key LzLi (or the tab key L“±lI ) to move to the next input field. 


d (n) 


d i r-j \ 



|r-:M 


The second input field is for the order of the derivative. If left empty, a first order derivative will be 
computed. This is what we want, hence we need not input the number one, but can move to the next 

input field with S or 


fit" 


There are several techniques to input expression x 2 ln(x) as the third argument. Note that it was en- 
tered already in the previous step. Hence, we could type ans (which is an abbreviation of “answer”) 
as a reference to the previous answer, i.e. the desired expression. As an alternative we use the up- 

arrow key m 

to highlight the last answer, . . . 

1 g| 



. . . then we use the enter key to paste a copy of the highlighted expression at the position where 

the input cursor was before choosing this expression. Another use of the enter key 1 completes the 
input. 



4(* 2 'lnW) 

2 , x , ln(x)+x 

dx 
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OS 

In general Derive ’s CAS engine is still more powerful than TI-Nspire ’s CAS engine, but there are 
already a few exceptions, such as solving trig equations, . . . 


so 1 ve (2 ■ cos (x)+ (s in (x) ) ^ = 1 j x) x = ^ ^ K 


... for which TI-Nspire may offer all solutions (by introducing a parameter, n2 in our example). 


OS 

Recently there has been a debate about how Derive approximates expressions involving inverse trig 
functions in degree mode. When we use Derive to simplify, everything is fine: 



However, when we approximate these expressions, the resulting angles are given in radians - even if 
the system is set to degree mode. 



One can obtain degree answers when using ARCSIN instead of ASIN, . . . 



r 1 1 


#12: ARCSIN 

— 



L 2 J 


#13: 


30 


r 1 1 


#14: ARCSIN 

— 



L 3 J 


#15: 


19.47122063 


. . . but this is not satisfactory when solving trig equations. 


#1: 

SINOO 

= 



2 

#2: 

Angle : 

= Degree 



r 1 ^ 

#3: 

NSOLVE 

SIN 00 = , x 

L 2 J 

#4: 




k = -3 , 665191456 
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If the above is a nuisance to you, you probably will find the answers provided by TI-Nspire to your 
liking. 



TI-Nspire offers the option to represent complex numbers in either rectangular or polar form: 


c Solve 


(v=i J x) 


-1 Js . -1 J 3 . 

x- — +- — i or x= — 1 or x=l 

2 2 2 2 


cSolve(.x^ = l J .x) 

X— ( 1 120)orx=(l I— -120)orx=l 


1 

TI-Nspire also offers a statistics menu with a lot of useful functions. 



One-Variable Statistics 

Two-Variable Statistics 

StatMatrix 

Regressions 

Random 

Distributions 

Tests 

Confidence Intervals 


03 

The user interface allows the typical Windows-standard mouse-click-and-drag method to select any 
part of an expression on the screen. 



In combination with the usual cut-copy-paste technique this is useful when building new expressions 
from old ones. However, other than in Derive, the selection of portions of an expression is not re- 
stricted to valid subexpressions. 


-1 

X— — + 
2 

a . 

•iorx= 

2 

-1 f . 

— 1 or x=l 

2 2 


5 r- 

— i or x= 

.1 1 
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Graphs & Geometry 

We use the lnsert>Graphs&Geometry command to add a page for function graphing and interactive 
dynamic geometry. 


Ei0O00i@ 


1)3 


|au (*■)— 

1 947 1 2 n 

cEcJwj (jr i - 1 

H 1 1- 

*-(l £. -I2D| or -1 

CEd^i (k 

1 2 

-i I 

2 2 1 

% 


DEC AJTQ tser 

11/39 



e,’ 



1 

1 3 


P 

n 

K 


1 ^ filth. 


The preview pane on the left now shows a reduced calculator page and a reduced graphing page, both 
of which are not interactive. One can easily flip between pages (i.e. making them interactive) by click- 
ing the desired page in the preview pane. 

In the upper right comer of the application window, next to the application window’s icon menu, the 
string 1)3 indicates that this is Page 3 of Problem 1. 

By the way, TI-Nspire offers a Presentation View and a TI-Nspire CAS+ Handheld View. The first one is 
very useful if one presents in a bigger room, the second one is very useful when the students use 
handhelds. In both views the preview pane is not visible. 

Insert Tools Help 


| | Layout ► 

TI-Nspire CAS+ Handheld View 


n 


Presentation View 


Shows an enlarged version of the page 


Back to function graphing. To obtain the graph of y = x 2 , we enter the function term x 2 at the bot- 
tom of the graphing window, hence turning the system reserved function name fl into a user-defined 
function. As it is customary in TI’s graphing calculators,//, are reserved for this purpose. This 

technique is similar to producing graphs on a graphing calculator, where one enters the function terms 
in a “y=”-window. In TI-Nspire one sees only one function (and its definition) underneath the graph- 
ing window and one can use the up CD and down [T] 

arrow keys to scroll. 


Ejy /j(x)=*-^ 


Also here we have the two-dimensional input of expressions. After entering the input with L_ _ J, the 
graph is drawn and is appropriately labelled with f1(x)=x 2 . 
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This graphics easily can be changed in many ways simply with the mouse. We choose the Pointer tool 
from the first menu, grab the horizontal axis at its left end (grabbing is performed by moving the 
mouse to the desired position, then clicking and holding the left mouse button) and extend it further to 
the left, . . . 


V 

'ip--, 

7 

-—1 

-6 

1 


. . . then we grab the horizontal axes at one of its tick marks and change this axis’ scale factor, . . . 


* V 

ly 

-4.5 

0.5 


. . . and finally we grab the entire plane (it can be grabbed on any white spot) and shift it around the 
window. 



7 

-3 

0.5 


For certain classes of expressions one can grab the graph, change its shape, and observe how the func- 
tion term changes accordingly. 


\ : 


\o5 

/ fl(x)=0.59'(x) 2 

-3 

0.5 



Using the Attribute command, one can change the appearance of a graph or how the function graph is 
displayed. 



4.5’ 

1 

\ 

\ 

V 

\ 

\ 

ft 

i 

1 

/ 

i 

i 

M 

# 

# 

t 

f \2 

(3/3) Line style is dashed / ^ 

\o.s 

s. ^ 

/ T1 y 

_ / 

ro=. 


-3 

0.5 

V 
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TI-Nspire also offers many of the dynamic geometry features known from Cabri Geometre such as 
drawing a circle: 


V 

\o,5 

s ^ 

y /fl(x)=0.61'(x) 2 

-3 

“G> 


Lists & Spreadsheet 


We use the lnsert>Lists&Spreadsheet command to add a page connected to a spreadsheet. This spread- 
sheet offers both column logic and cell logic. One can do the standards of what one can do with stan- 
dard numerical spreadsheets; in addition this spreadsheet can carry and process algebraic expressions. 
The only limitation may be that within the cells the expressions are displayed in linear format, only by 
selecting a cell one can view its contents in two-dimensional format at the bottom of the page window. 



This is ideal for having students recognize patterns. Once they have an idea, such as for the above, the 
teacher may challenge them by adding an example for which the validity of their pattern needs some 
additional thinking, such as in the following screen image. 



x A 20 

20*x A 1 9 

x A 21/21 



6 

|x A (-3) 

-3/x A 4 

-1/(2*x A 2) 




Notes 

We use the lnsert>Notes command to add a page for producing text. 

This text processing application offers 
standard features such as bold, underlined , 

sub scri P tedj SUperscri P te d etc, 

It allows for many special characters 
such as for denoting angles, Z_(ABC), 
triangles, A(ABC), segments, AB, and more. 

If one refers to definitions made in 
other applications, such as f 1 (x), these can 
be selected and evaluated: f1(x)=J 
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By applying the Evaluate selection command, the selected expression is algebraically simplified, i.e. 
also the text processing application is connected with the computer algebra application. 


If one refers to definitions made in 
other applications, such as fl(x) # these can 
be selected and evaluated: f1(x)=x A (2)| 


Miscellaneous 

Within a TI-Nspire document, user-defined functions and variable assignments are local to a “Prob- 
lem”, within which they are valid in all pages. Other problems in the same document have their own 
local variables. 

A very useful command is the general Undo command, which undoes whatever has been done last. A 
corresponding Redo command is also available. 

tfj pi 


For fast access to certain symbols or build-in functions, one can switch on the Keypad, which can be 
put next to the application window for convenient input. This is how the calculator keyboard looks 
like! 



Other sources for supporting the inputting of expressions are the Catalog with a list of all build-in 
functions in alphabetical order including a quick “online help” in form of “expressive” templates at the 
bottom of the window, . . . 



. . . and a Symbol Palette for special symbols and input templates. 
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Closing Remark 

TI-Nspire’s innovative double-platform design (handheld + 100 % compatible computer software) and 
the seamless integration of the various components (CAS, graphing, geometry, spreadsheet, text proc- 
essing) make it a very powerful teaching and learning aid for teachers and students that is unique on 
the market and has great potential. Please help testing it. I’ll be happy to collect the feedback and for- 
ward it to the appropriate people. 


Beyond Polynomial Regression and 
Beyond quart reg on the TIs! 

Josef Bohm, Wiirmla, Austria 

Well known - since several years even in school mathematics - are linear and quadratic regression, 
both based on the least squares method. One of the reasons that regression models found their way in 
class room is the increasing importance of statistics and another one is the availability of tools which 
make the tiresome numerical calculation easier or even unnecessary. 

Given is a set of bivariate data {(xi ,yi), (*2,72), > (x n ,y n )}. Presentation of the data in a scatter dia- 

gram sometimes leads to the conjecture that there might be a linear relationship between x- and 
y- values. Regression tries to find a linear model function/(*) = kx + d such that the sum of the squares 
of the differences between real and theoretical function values becomes minimal: 

n n 

= X(/( x z) — k/) 2 = X(^ x z + d - y t ) 2 = Minimum 

z = 1 z = 1 

This is an extremal value problem with two variables k and d. The values of them can be found by 
partial differention of £l(k,d) with respect to k and d and then solving the resulting system of linear 
equations for k and d. 

This method is going to be generalized. Let's have a demonstration example: 

(I will start working on the handheld, the DERIVE realisation will follow.) 


Xi 

1,36 

2,14 

2,76 

4,40 

4,99 

5,85 

6,89 

8,25 

8,50 

10,09 

y i 

3,65 

2,65 

-0,20 

-2,12 

-1,96 

0,49 

3,44 

2,88 

0,97 

-2,11 


Store the data in two lists X I und y I and have a first impression by inspecting the scatter diagram: 
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We can be sure that linear regression is not appropriate in this case. We could try polynomial function 

of higher degree, but the scatter diagram reminds us on trig 
functions, doesn't it? 

So we look for an approximating function 0(x) defined as 
a linear combination of the set (list) of functions: 

{1, sin(x), cos(x)}. 

O(x) = a • 1 + b • sin(x) + c • cos(x) 

We obtain an extremal value problem with three variables a , b and c: 

®( x i ) = *1 f\(xi) + a 2 / 2 (Xi ) + a 3 / 3 (x-) 

N , : A v 

Q(a,9,c) = a - 1 + &*sin(x.) + c* cos(x / ) - y.) = Minimum 

i = 1 

9Q ^ 

= ^(2 (a + b sin(x z .) + c cos(x z ) - y.) • 1) = 0 

da i = i 

0Q ^ 

= X (2 ( a + b sin(x . ) + c cos(x z . ) - y t ) • sin(x z . )) = 0 

db i - 1 

9Q N 

= X (2 ( a + b sin(x z . ) + c cos(x z . ) - y . ) • cos(x z . )) = 0 

9c i = i 

The resulting equations can be cancelled by 2, expanded and rearranged: 

a y \ • 1 + b y'sinx i • 1 + c^ j cosx i • 1 = *1 

N N N N 

a y\ • sin x. + b'^ j sinx i - sin x t + c^cos x t • sinx t = ^y. • smx. 

N N N N 

a yi • cosx i + b^ j sinx i - cosx i + c^cos x f - cos x t = ^y i - cos x t 

N N N N 

The system is linear again. Our task is to set up the system for the given data (7V = 10) and solve it. 
Then add the graph of the “regression line” to the scatter plot. We can use the lists X I and y I . 

Note: S U m( I i S t ) gives the sum of all elements of a list. 

10 a + 1.533626 b - 1.336507 c = 7.690000 

1.533626 a + 6.006008 9 - 0.829608 c = 16.119240 

-1.336507 a - 0.829608 b + 3.993992 c = 2.871270 

which leads to: 0(x) = 0.544245 + 2.748197 sinx + 1.471855 cos x. 




The Calc-Menu £ in the DATA-MATRIX-Editor of the 
TIs offers option Si n R e g which gives a sinuoidal regres- 
sion line. 

The result is: 

y(x) = 3,09219 sin(l, 02073 x + 0,37917) + 0,54421 
thick: 0(x), line: y(x) and dot: quart r eg 
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We compare the quality of both regression lines by calcu- 
lating the sum of the squares of the residuals. The Tl-built- 
in trig regression proves to be a little bit better. 

Calculating the trig regression line is much much more 
complicated (see DNL# 46) 

The method shown above can be performed very elegant by using matrix algebra. First of all we create 
matrix M, the so called Normalmatrix: 


| t f — |fl 1 gebrajCaic |OiherJPr gm 1 0 |C 1 ean Up I I 

U - 1 i 336507 - . 829608 3.993992 

1 . 000000 0. 000000 0.000000 .544245 ' 

0.000000 1.000000 0.000000 2.748137 

Lo. 000000 0.000000 1.000000 1.471855. 

■ £un((y6CxlJ - yl) 2 ) 1.772076 

■ sun((y7Cxl) - yl) 2 ) 1.644764 

sumC C»j7<x 1 >-«jl > A 2 > 


'/l(*l) fl( X 2 ) fl( X i) /l(*lo) > | 1 1 1 1 

M = / 2 (x,) f 2 (x 2 ) / 2 (x 3 ) f 2 (x jo ) = sin(xj ) sin(x 2 ) sin(x 3 ) sin(x 10 ) 

v / 3 Oi) M x i) M x 3) / 3 ( x io)J ^cos(xj) cos(x 2 ) cos(x 3 ) cos(x 10 ) y 

What ist he result of M -M‘? 

Solving the simultaneous equation is easily done: 

Please compare with the solution from above! 



Her is a problem for you(r students): 



The approximating function shall be a linear combination of: 

f(x) = 1 ,/ 2 (x) = x 2 ,/ 3 (x) = (x+1) 1 und/ 4 (x) = (x+2)' 2 . 

Try to use appropriate Tl-functions to keep edit- 
ing as little as possible. 

Y ou can read off the solution: 


Possible questions for students: 

Can we use a polynomial of higher degree as regression line? Do you find another sensible 
set of basic functions? Compare the quality of the regression models by calculating the SSE 
(sum of squared errors = sum of squares of residuals). This method can help to overcome 
the restriction of polynomial regressions on the Tl. Do you have any idea how to do? 
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Realisation with a DERIVE-function 


Transfer to DERIVE is very easy, because we can use the TABLE-function very efficiently to create 
the normalmatrix. A three-lines program does the whole job: 


genreg(fl ist , pts, nm, rs) = 

Prog 

#1: nm = REST (TABLE (fl ist , x, ptsul) 1 ) 

rs = [VECTOR (ptsi i2*nm|k, k, D3M(f! ist))] 1 
((nm*nm 1 ) A (-l)*rs)ul*flist 


' 1.36 

3.65 ‘ 

2.14 

2.65 

2.76 

-0.2 

4.4 

-2.12 

4.99 

-1.96 

5.85 

0.49 

6.89 

3.44 

8.25 

2.88 

8.5 

0.97 

. 10.09 

-2.11 . 


Big surprise? Or not? 

FIT([x, a + b*SIN(x) + c*COS(x)], test) 


|1, 47185 -COS(x) + 2 . 74818 »SIN(x) + 0.544245| 


Please read Karsten Schmidt's com- 
ment on this kind of regression (page 5) 
and read the paragraph on the FIT- 
function in the Online Help of 
DERIVE. 

You might ask if our work was super- 
fluous and waste of time? My opinion 
is: No, because we could make a “Black 
Box” to a “White Box” with relatively 
little effort. 


#3: genreg([l, SIN(x) , COS(x)], test) 

#4: 1 . 471854558 *C05(x) + 2 . 74S1S7037 *SIN(x) + 0.5442452255 


Let's compare with the si nreg(datamatri x) program for trigonometric regression. 

In 2002 Don Phillips provided a wonderful tool based on the Levenberg-Marquardt and Gauss- 
Newton Methods for nonlinear regression models. 3 


/ 5: 

# 6 : 

*7: 

*3: 


sinreg(test) 

' Regressions! in-ie: y = 0.544210 - 3 . 09218. SIN(l . 02073 *x + 3.52075) 

SummeRes A 2 : 1.64476 

TABLE(1. 47185. COS(x) + 2. 74818. SIN(x) + 0.544245, x, -1, 12.5, 0.01) 
TABLE(0. 54421 - 3 . 09218. SIN(1 . 02073 .x + 3.52075), x, -1, 12.5, 0.01) 



3 Derive Newsletter #46, Statistics-Tools for DERIVE and the TI-92 



D 


-N-L#63 


Josef Bohm: Beyond Polynomial Regression 


p29 


test 2 is the data matrix of the 2 nd example (and try again with the FIT-function and compare!): 



- F jT 

2 1 1 


EXPAND 

genreg 

1, x , 

, test2 



* + 1 2 




L U + 2) J 



94.79609995 27 . 22541J55 2 

+ + 0.06679516555 x - 15.948999J4 

#7: 2 x + 1 

(x + 2) 



Programming for the CAS-TI-Calculator 

genreg (f uncti onl i st , data matrix) is the respective program fort he TIs. 

A special list of functions (powers of x) leads immediately to a higher degree polynomial regression 
pol reg (degree , data matrix). (See the screen shot on the next page.) 

genreg 
(fl i st , pts) 

Prgm 

Local nm,rs,xv,yv 

pts T [1 ]^xv: pts T [2]^yv 

newMat (dim(f 1 i st) [2] , dim(xv) [2] )->nm 

newMat (dim(f 1 i st) [2] ,1 )->rs 

For k, 1 ,dim(fl ist) [2] 

For 1,1 ,dim(xv) [2] 

1 imi t (f 1 i st [1 , k] , x , xv [1 , 1 ] )->nm[k, 1 ] 

EndFor 

dotP(yv , nm[k] )->rs [k, 1 ] 

EndFor 

(nm*nm T ) A ( "1 ) *rs->rs 
dotP(rs T [1 ] , f 1 i st [1 ] )->regf 
Di sp regf 
Disp "" 

Disp " Gespei chert unter regf" 

EndPrgm 


I'FiTTif f£t y f::t y fht y fs Y FfiT T i 
|fllgebra|Calc|0ther|PrgmI0|ciean Up| | 

■ 

T736 3.65 
2.14 2.65 

2.76 -.2 

4.4 "2.12 

4.99 "1.96 

5. 85 .49 

6.89 5.44 

4 data 

1 . 3600 3.6500 
2.1400 2.6500 

2.7600 -.2000 

4.4000 "2.1200 

4.9900 "1.9600 

5.8500 .4900 

6.8900 3.4400 


...nresCEl .sin<x> .cos(x)] ,data>l 

MAIN FiliD HFFRDK FUhK 50/50 
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4 1 

- 94.796100 27.225414 

068795 x 2 

(x + 2 . 000 GG 0) 2 x+ 1.000000 

Stored as regf 


MAIN RAD APPRUK FUNC £9/30 1 



Compare the results! 


The Tl-program for polynomial regression pol 
reg ( ) refers to genreg ( ) : 




F£r Yf::tTfHtT FF Y FfiT Y 
Control I/O Oar Find... Mode 


: po 1 reg (. g , pts > 

: Frgn 
: Local fl 

: augnentC £1} , seqtx^k, k ? l ? g>>-^fl 
: lisUnaUflMl 1+fl 
: genregtf 1 ? pts) 

: EndPrgn 


£3 


Y F]:t T FHt Y fl Y" Ffi' 

fllgebra|Calc|Other|FrgnIO|Clean 


i 


Done 


■ regf 

27 225414 r> 

i > + x + 1 . 000000 + ■ 066795 ■ x2 - 15.948999 

■ polreg £6 , data 2 ) Done 

■ regf 

< 90462241 x 2 - 3 . 935 Q 2253 -X - 6.28012454 

resrfl 

MAIN RAD AF'F'RDH FUNC £5/30 


MAIN 

RAD APPRUK 

FUNC 




ss 

. f; -r V T-J'^ 

LVs^L;r-i:|t. 5 i C | : .;th(-E 

'IprgnlO^L 



.00000151 x 6 - . 00014535 x 5 + . 00530918 x 1 


Stored as regf 


MAIN RAD AFFRDK FUNC £H/30 


Compare again with the DERIVE result applying the FIT-function: 

r 6 5 4 3 2 1 

#9: FITCLx, a-x + b-x + c-x + d-x + e-x + f-x + gj, test2) 

6 5 4 3 2 

#10: 0.0000015 OS 06 -x - 0.000145352.x + 0.00530917-x - 0.0931230-x + 0.904622-x - 3.93302-x - 6 . 2S012 


A Closing Example: 

A Lorenz Curve is defined by the points (0|0), (0.2|0.03), (0.4|0.1), (0.6|0.25), (0.8|0.55) und (1 1 1) 
(= Quintil-values). The curve shall be modelled by a quartic passing the origin. 


This is no problem in DERIVE: 


#17 : lorenz ;= 


/IS: FIT (Lx, a-x 


' 0 

0 

0.2 

0.03 

0.4 

0.1 

0.6 

0.25 

O.S 

0.55 

. 1 

I . 

4 

- x + 

3 

b-x + 


4 3 2 

/19: - 0. 5156457503 -x + 1 . 00S615537 -x - 0.6003001062 - x + 0.207S1S7250 x 
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1 Jf 

o.s j ■ 

-0.6 / 

0.4 / 


0.2 0.4 0.6 O.S 1 
-0.2 

• Note: In economy the Lorenz curve is used to represent the inequality of eg. income distribution 
among different groups in the population of a country. This inequality of the distribution is meas- 
ured by the GINI coefficient, which measures how far the Lorenz curve falls below y = x. It is the 
ratio of the area between y = x and the Lorenz curve and the area between y = x and y = 0. (This is 
a very welcome application of the defined integral.) 


From today on this is no problem for you as a Tl-User, too: 


rriT^T fzt y fjt y fh^ ^ 
| t |fl 1 gebra |Ca 1 c lother 

r fe Y FfiT T 1 
PrgmIOlClean Up| 




L 1 . 00000000 1 . 

00000000J 


'0. 

0. ' 


'0.000000 

0.000000] 


.2 

.03 



200000 

. 030000 


.4 

. 1 

+ lorenz 


400000 

. 100000 

■ 

.6 

.25 


600000 

. 250000 


.3 

.55 



S00000 

. 550000 


.1. 

1. . 


.1 

. . 000000 

1 . OOOOOOJ 


- .25] [ 

.8 . .55] [1 


,1 - ] ]-KLorens 

MAIN 


rad apprdk 


FUNC £9/30 



iTupf~ 


I .3 .551 

I .500000 

. 550000 

[l- 1- J 

L 1 . 000000 

1 . 000000J 

■ genreg([ x x * 1 2 x 3 x 

■ regf 

4 ] ? lorenz) 

Done 

-. 515646 x 4 5 6 7 + 1. 903616 x 3 - . 600390 x 2 > 

■ ". 5156457503 x 4 + 1. 

9086155382- 

x 3 - .60* 
Done 

ans<l>-*y3<x>l 

MAIN RAD APPRDX 

FUNC 30/30 




For higher degree polynomial regression it can happen that the matrix will become singular caused by 
internal calculation inaccuracy. 
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Contour Plots - Level Curves 


Peter Schofield, UK & Josef Bohm, Austria 

In DNL#59 appeared a request from Manuel Garcia Mateos for plotting level curves 
(= contour plots) with DERIVE. Manuel was right complaining that other major CAS- 
programs (MAPLE, MATHEMATICA) provide packages for creating contour plots either 
in the xy-plane or as a 3D-plot. Producing a plot in the xy-plane is not so difficult using the 
VECTOR-command and DERIVE 's ability to create implicit plots. As a “Graphic Freak” I 
wanted to overcome the deficiency of presenting the contour curves in 3D and started 
thinking about a program. Now, since I am able to use three thicknesses of lines in the 
plane we have more varieties plotting even in plane. Peter Scofield who is also a graphics 
enthusiast and a real graphics wizard joined the party and gave his version of facing the 
challenge. Moreover Peter produced a toolbox for plotting implicit forms of three variables 
and found an improved version of NSOLUTIONS and NSOLVE. 

I'll demonstrate what we did using function f: 

Finding the contour plots in xj-plane is not too difficult: 

2 2 

#1: 2 2 1 - x - y 

f := (x + 3*y )*e 

3 


0 


Two more examples of plane contour plots with their graphs on the next page. 

862 5442 3 226 

#15: f2 := 7(16. x + 8»x .(8. y + 17) + 112. x + x .(96.y + 136. y + 801) + 28. x .(17 - 8.y ) + 2»x .(32. y 

42 42 8642 

- 68. y - 1247. y + 1186) - 28«x.(12.y - 17. y - 64) + 16. y - 136. y + 801. y - 1980. y + 4096) 

#16: VECTOR (f 2 = z, z, 0, 200, 20) 




-y 4*y 2*y 2 

c -7(c + 2*c *(1 - 2*COS(x) ) + 1) 

#23: f3 := 


2 

#24: VECTOR (f 3 = z, z, 0, 1, 0.1) 

We want to see the level curves in 3D-plot - either on the surface or as projection on the xj-plane - 
and we would like to plot the 2D-level curves in various thicknesses. 



For the following I need my auxiliary function mysol uti ons, which - as I am hoping - finds all nu- 
meric solutions of equation equ vor variable equ_v in the range between x = start and x = end. 

#3: mysol ut ions (equ, equ_v, start, end, step := 0.1) := APPE N) ( VECTOR (NSOLUT IONS (equ, equ_v, a_, a_ + step), 

a_, start, end - step, step)) 


This is the first half part of my programm con tou r () . 

contour (f, level, xst, xend, yst, yend, xstep := 0.2, ystep := 0.2, acc := 0.2, equ, x_, y_, pts) := 
Prog 

pts := [] 
x_ := xst 

f := SUBST (LHS(f) - RHS(f) , z, level) 

Loop 

equ := SUBST (f, x, x_) 

If equ 

pts := APPEM)(pts , [[]]) 
pts := APPEMD(pts , [[]]) 

pts := APPEhD(pts, VECTOR ([[x_, v, level]], v, mysol ut ions (equ , y, yst, yend, acc))) 

/S: x_ :+ xstep 

If x_ > xend exit 
y_ := yst 


I enter any equation with three variables x, y and z as f . 1 evel is the distance of the level plane from 
the xy-plane. The plot range is defined by xst < x < xend and yst <y < yend. In the 3 rd line f will 
change to an implicit form f(x,y) = level. In the following loop we substitute x by all values from xst 
to xend with an increment of xstep (default value 0.2) and try to find all solutions fory, which leads 
to points on the surface - if there are solutions. The first two parts in the IF-construct are necessary for 
special cases (f from line 3 can be returned as f al se, . . .). If you are fine with the default values, then 
entering parameters ends with yend. acc sets the accuracy for the mysol uti ons-procedure. 

Then we repeat the same procedure for the sequence of y- values to have a tighter point set. 

See my first results (don't mind some calculation time!) 

*9: contour(f = z, 0.25, -3, 3, -3, 3, 0.1, 0.1, 0.5) 

#10: VECTOR (contour (f = z, z, -3, 3, -3, 3, 0.1, 0.1, 0.5), z, 0.5, 2.5, 0.5) 

Highlight the expressions, switch to the 3D-Plot Window (set Options > Approximate before Plotting) 
and Insert > Plot. When DERIVE is ready (in the background) you will be asked for the Point Size 
(Small / Medium / Large). Try any of them. Then you will be asked for the Plot Colour. Set Scheme 
Custom and set Top Maximum and Minimum the same colors. 
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Contour plots of #15 from above: 




This is not so bad. But if you choose Point Size Small, the points are not really small. I wanted to have 
only dots as points. It is not possible to have discrete points appearing as dots - either the points are 
connected or if dicrete, then small (which is not really small), medium or large. (This is what I know). 
So I wanted to have a “contour dust”. Peter gave the idea: plot a segment of length 0. [1,2, 3; 1,2, 3] 
shows a tiny dot. So I only had to change the output form of contour to get contourdotO with the 
same parameters as before. 

#12: VECTOR (con tourdot (f = z, z, -3, 3, -3, 3, 0.1, 0.1, 0.5), z, 0.25, 3, 0.25) 

*22: VECTOR (con tourdot(f 2 = z, z, -3, 3, -3, 3, 0.1, 0.1, 0.5), z, 0, 200, 20) 

5 5 5 

*36: VECTOR ( |x | + |y| + |z| = 32, z, 0, 2, 0.1) 

5 5 5 

*37: VECTOR(contourdot(|x| + |y| + |z| = 32, z, -2, 2, -2, 2, 0.1, 0.1, 0.5), z, -2, 2, 0.2) 


See the plots on the next page: 
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And finally there is contourdotpl ane() which gives the projections of the level curves in the 
xy-plane but as a 3D-object: 

#14: VECTOR (contourdotpl ane (f = z, z, -3, 3, -3, 3, 0.1, 0.1, 0.5), z, 0.5, 3, 0.5) 



As it was mentioned in the contribution about plotting thick lines, it was not too difficult to adapt the 
basic idea creating contour plots for doing implicit plots with thick lines (which are in fact only thick 
points in a row). 

impplotO is a adaptation of contourdotplane(). 

Peter sent a mail together with his files pol ycontou r . mth and pol ycontou r . dfw and just recently an 
improved version: 

#1: Procedures (based on Josef's contour idea) for 3D-plotting contours of surfaces of form f(x,y,z) (=0) 

or f(x,y,z) = g(x,y,z) using a relatively quick N_SOLLTTIONS procedure. 

N_SOUJTIONS ( f , x, a, b, c, i, z, p := []) := 

Prog 

Peter found a quite better auxiliary program for finding the numeric solution of an equation. I'll pre- 
sent this programs N_S0LUTI0NS() and N_S0LVE() in an extra paragraph. 

#3 : 2D-CONTOUR-POINT-PLOTS 

ContourPts_2D(f , z_, stx := -4, endx := 4, sty := -4, endy := 4, stepx := 0.1, stepy := 0.1, s, f_, p := []) := 
Prog 

Approximate and 2D-Plot 

2 2 2 

#5: VECTOR (Con tour Pts_2D(x + y + z - 16, 1 ) , 1 , 0 , 4 , 0 . 5) 

2 2 2 2 

#6: VECTOR (Con tour Pts_2D((x + 3*y )*EXP(1 - x - y ) - z, z_, -2, 2, -2, 2), z_, 0, 3.2, 0.2) 
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For 3D-representations Peter offers ContourPts_XY() and ContourDots_XY(). Using the first pro- 
gram you can choose between the three point sizes Small / Medium / Large and the second one gives 
us the “dust”. An important EXTRA is the fact that his function can process Boolean expressions and 
operators. 


#8 : 3D-CONTOUR-POINTT-PLOTS 

ContourPts_XY(f , z_, stx := -5, endx := 5, sty := -5, endy := 5, stepx := 0.1, stepy := 0.1, s, f_, p := []) := 
Prog 

#10: ContourDots_XY(f , z_, stx := -5, endx := 5, sty := -5, endy := 5, stepx := 0.1, stepy := 0.1) := ADJOIN([], 
VECTOR([u, u], u, REST(ContourPts_XY(f , z_, stx, endx, sty, endy, stepx, stepy)))) 


#11: [Precision := Approximate, PrecisionDigits := 4, Notation := Decimal] 

#12: Approx, and 3D- PLOT : 

2 2 2 2 

#13: VECT OR (Co n to u r Pt s_X Y ( z = (x + 3.y ).EXP(1 - x - y ), z_, -3, 3, -3, 3), z_, 0, 3.2, 0.2) 

#14: It can also process Boolean expressions and operators!! Approx, and 3D-PLOT : 

2 2 2 

#15: VECT OR (Co n to u r Pt s_XY (x + y + z = 16 A x.y > 0, z_, -4, 4, -4, 4), z_, -4, 4, 0.2) 




x 
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In my short introduction I announced that Peter 
provided a file for doing implicit 2D and 3D-plots. 
Here are some mouth watering examples. I'll 
come back to that in the next newsletter. 

Among the files which you can download your 
will find contours J.dfw, contoursj.mth, polycon- 
tours. dfw and polycontours. mth. The dfw-files 
contain the worked examples and the mth-files 
are containing the functions only. They should be 
loaded as utility files in the background. 

Happy scanning, 

Peter & Josef 

* 66 : 



5 5 5 5 

Lnp'Hc'itDots(|x| + |y| + |z| - 2 , -2, 2) 



mysol uti ons versus N_SOLUTIONS 


Josef's version for numeric solutions of equation equ (first lines of code): 

mysol ut ions (equ, equ_v, start, end, step := 0.1, list, sol_list) := 

Prog 

list := APPE hD ( VECTOR (NSOLUT IONS (equ, equ_v, a_, a_ + step), a_, start, end - step, step)) 
If Dm(list) = 0 
RETURN list 

sol_list := [FIRST(1 ist)] 

/I : Loop 

Peter's version for numeric solutions of equation f (first lines of code): 

N_SOLirriONS(f, x, a, b, c, i, z, p := []) := 

Prog 

C := [[a, b]] 

Loop 

■i := DUI(c) 

Loop 

a := (FIRST 
b := (FIRST(c))i2 
c := REST (c ) 

If a < b 

*2 : Prog 

z := NSOLUTIONS(f, x, a, b) 

There is also a N_SOLVE-function which might substitute the DERIVE NSOLVE. Inspect the file. 
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See first, what DERIVE is offering:: 

( - 0.2.x 


#3: 

NSOLUTIONS|^5 »e 

#4: 

[-5.65258199] 

#5: 

SOLUTIONS^, c 

#6: 

[-9.424133421, 


•SIN(5.x) 


•SIN(S.x) 


5 T' X ) 


#6: [-9.424133421, «, -<», + 47 . 1238898. L , « + 15 . 70796326. L , « - 15 . 70796326. L , « + 31 .41592653 , « 

- 31. 41592653. l] 

If you add ‘Real’ as last parameter, then we obtain -9.4241 ... as the only one solution . . . 

Now let's try what we can do: 


:ions^5 


- 0.2.x 1 ^ 

#7: mysolutions| 5 *c *SIN(5*x) = , x, -10, lol 


*8: [-9.424133421, -8. 79724213, -S . 167184593 , -7.540996263, -6.910050444, -6.234996034, -5.65253199, 

-5.02945696, -4.39444973, -3.774392193, -3.134774639, -2.522350526, -1.370223676, -1.230367643, 
-0.5633236745, -0.2167233973, 0.2296115766, 0.5433226014, 1.296334646, 1.35356236, 2.539524499, 
3.11760029, 3.792477949, 4.37625373, 5.043333594, 5.632915039, 6.305621164, 6.333424393, 
7.563335603, 3.14303332, 3.323010163, 9.396307952] 


31 sec 

f - 0.2.x 1 ^ 

*11: N_SOLlTTIONS 5*c .SIN(5.x) = , x, -10, 10 

#12: [-9.424133421, -8.797242130, -8 . 167184593 , -7.540996268, -6.910050444, -6.284996034, -5.652581990, 

-5.029456960, -4.394449780, -3.774892193, -3.134774689, -2.522850526, -1 . 870228676, -1 . 280867648 , 
-0.5638286745, -0.2167283973, 0.2296115766, 0.5438226014, 1.296884646, 1.853562360, 2.539524499, 
3.117600290, 3.792477949, 4.376253730, 5.048338594, 5.632915089, 6.305621164, 6.888424393, 
7.563885608, 8.143038820, 8.823010163, 9.396807952] 


7.4 sec 


#13: 


sols := N_SOL UT IONS 5»c 


0.2.x 

»SIN(5 »x) 


1 

, x, -10, 10 

x 



I took some other equa- 
tions and finally I must 
surrender and admit that 
Peter's tool is much better 
than mine, congratula- 
tions. 

The nice thing is that 
Manuel Garcia's request 
on contour plots lead us to 
fine results. We hope that 
some of you will find them 
useful for their work. 
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Find intersection points of a parametric 3D ellipsoid and a 
3D parametric straight line 

Bill Wilburn, USA 


solutions=N [Solve [x==4 Cos[s] Sin[t]&&y==3 Sin[s] Sin[t]&&z==2 Cos[t] 
&&x==u+l&&y==3 u&&z==u+l ,{x,y,z,s,t,u}]] 

Solve: :ifun: Inverse functions are being used by Solve, so some solutions 
may not be found; use Reduce for complete solution information. More... 

{ { x^O . , y^-3 . , z-»0 . , u— >— 1 . , t->-l . 57 08 , s->l . 57 08 } , { x^O . , y^-3 . , z->0 . , 
u->-l. , t-»1.57 08,s->-1.57 08}, {x^l . 52 381, y->l . 57143, z^l . 52 381 , u^O . 52 381 , 
t->-0 .7 04547, s^-2 . 19959} , {x^l . 52 381, y->l . 57143, z^l . 52 381 , u^O . 52 381 , 
t->0 . 7 04 54 7 , s^O . 942 } } 

{x,y,z}/. solutions 

{ { 0 . , -3 . , 0 . }, { 0 . , -3 . , 0 . }, 

{1.5238095238095237,1.5714285714285714,1.5238095238095237}, 

{1.5238095238095237,1.5714285714285714,1.5238095238095237}} 

plotl=ParametricPlot3D [ {4Cos [s] Sin[t] ,3Sin[s] Sin[t] ,2Cos[t]}, 

{s , -7T,7r} , { t , -7T,7r} ,Axes^True , AxesLabel->{ "X" , M Y" , "Z" } ] ; 
plot2=ParametricPlot3D [ {u+1 , 3 u,u+l} , {u,-2,2}] ; 
plot3=Show [Graphics3D [ {RGBColor [1 , 0 , 0] , PointSize [0 . 025] , 

Point [ { 0 , -3 , 0 } ] } ] ] ; 

plot4=Show[Graphics3D[ {RGBColor [0,0,1] , PointSize [0 . 025] , 

Point [{1.52380952380952, 1.57142857142857, 1.52380952380952}] }] ] ; 

Show [ {plotl ,plot2 ,plot3 ,plot4 } ] 

Null 
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Parameter form: 

#1: [x_ := 4.C0S(s).SIN(t), y_ := 3 .SIN(s) .SIN(t) , z_ := 2*C0S(t)] 

Implicit form: 


*2: 


# 3 : 


2 2 2 

X_ y_ Z_ 

+ + 


16 


9 


#4: SOLUTIONS 


#5: 


x y z 

+ + =lAx=u+lAy=3*uAz=u+l, [x, y, z, u] 

16 9 4 


0-3 0 


-1 ' 

32 11 32 

11 

21 7 21 

21 . 

i — i 

o 

m 

o 

l_l 



' 32 11 32 ' 



. 21 7 ’ 21 . 



[0, 


I, 0] 


# 6 : 


#7: 

[1.523309523, 1.571423571, 1.523309523] 

#3: VECTOR([[u + 1, 3*u, u + 1]], u, -3, 3, 0.01) 

I use the parameter form for the 3D-Plot (no mesh lines) and a similar trick as in the contribu- 
tion on "Thick Lines" to plot the straight line as a thick line (Point Size Small is sufficient). 

The intersection points have Size Large. 



We can do it in another way, too: 

Substitute sin(t) = v and sin(s) = w and then square the components of line and surface. 

The system is now consisting of polynomials only and can be solved by DERIVE. 

#9: 5UBST([x_ := 4*C05(s) *SIN(t) , y_ := 3 *SIN(s) *SIN(t) , z_ : = 2*C05(t)] , [s, t], [ASIN(v) , ASIN(w)]) 

#10: [x_ := 4*w*^(l - v ), y_ := 3*v*w, z_ := 2*^(1 - w )] 

[ 2 2 2 2 2 21 
#11: SOLVE ( L(u +1) = x_ , (3 • u) = y_ , (u + 1) = z_ J , [u, v, w]) 


p42 


Bill Wilburn: Special CAS-Guest MA THEM A TICA 


D-N-L#63 


#12: |u = -1 A V I 1 A W = 1, U=-l A V = 1 A W = —1 , U=-lAV=-lAW=l, U = -1 A V = -1 A W = -1 , U 

11 11.^185 4185 11 11.^185 4185 11 

= A V = A W = , U = A V = A W = - , U = A V = - 

21 185 21 21 185 21 21 

11.7185 4185 11 11*4185 4185 ' 

A W = , U = A V = - A W = - 

185 21 21 185 21 . 

f \ 11 

#13: VECTOR [u + 1, 3 ♦ u , u + 1], u, -1, 

l L 21 

0 -3 0 

1.523809523 1.571428571 1.523809523 . 

But there is also a “Black Box” - way, which doesn't need more or less tricky manipulations. 
Students should know that there is a Newton Method for solving “uncomfortable equations” 
and they should learn that this method can be extended to more dimensions. In the DERIVE 
Online Help we can find via the index item ‘system of equations’ a reference to Equation- 
Solving. mth and there is a description of NEWTONS(u,x,xO) as a tool for approximative 
numerical solving of nonpolynomial systems. Let's try: 

#15: NEWTONS([u + 1 - 4 *CQS(s) *SIN(t) , 3*u - 3 *SIN(s) *SIN(t) , u + 1 - 2*CG5(t)], [u, s, t], [1, 1, 1]) 

111 
0.5643320325 0.9712704229 0.7125703502 

0.5240209059 0.9423926175 0.7044230923 

#16: 

0.5238095633 0.9420000063 0.7045473966 

0.5233095233 0.9420000403 0.7045474134 

. 0.5233095233 0.9420000403 0.7045474134 

We can observe that the method converges to a solution: u = 0.5238..., s = 0.9420... 
and t = 0.7045... 

To find another solution we have to change the initial values (which is now an initial vector 
for u, s, and t): 

#17: NEVTTONS([u + 1 - 4 .COS(s) .SIN(t) , 3*u - 3 .SIN(s) .SIN(t) , u + 1 - 2.COS(t)], [u, s, t], [-1, -1, -1]) 

leads to the same solution, but: 

#19: NEWTON5([u + 1 - 4 *CGS(s) *SIN(t) , 3*u - 3 *SIN(s) *SIN(t) , u + 1 - 2*COS(t)] s [u, s, t], [0, 0, 0]) 

needs about 14 iterations to converge for the second solution. 

-1 -4.71233393 4.71233393 

-1 -4.71233393 4.71233393 

-1 -4.71233393 4.71233393 

So we obtained the same intersection points as above. 

#21: VECTOR([u + 1, 3*u, u + 1], u, [0.5233095233, -1]) 

1 . 523309523 1 . 571423571 1 . 523309523 

0-3 0 

We could refine NEWTONS to make the output more friendly. In fact we are only interested 
in the last rows to see if and where the method is converging. 
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#26: 


APPEW)([[x, y, z, u, s, t]], (REVERSE(NEWTONS([x - (u + 1), x - 4 .COS(s) .SIN(t) , y - 3-u, y - 

3.SIN(s).SIN(t), z - (u + 1), z - 2.C0S(t)], [x, y, z, u, s, t], [1, 1, 1, 1, 1, 1]))) ) 

[1, 2, 3] 

x y z u s t 

1 . 523809523 1.571428571 1.523809523 0.5238095238 0.9420000403 0.7045474184 

1.523809523 1.571428571 1.523809523 0.5238095238 0.9420000403 0.7045474184 

1.523809563 1.571428691 1.523809563 0.5238095638 0.9420000068 0.7045473966 

We don not only see the convergence, but also all solutions - 6 equations and 6 unknowns. 
For the other solution you have to take [0,0, 0,0, 0,0] as initial vector. 


#27: 


I sent this treatment to Wilburn and a very few days after I received another mail with an im- 
proved ( = more complicated) version of his problem. The centre of the surface moved out of 
the origin: 


Clear [sol] 

sol=NSolve [x==l . 414213562 Cos[s] Sin [t] -1 . 719669914 Sin[s] 

Sin[t]+1. 478397839 Cos [t] &&y==2 . 449489742 Cos[s] Sin [t] +2 . 217596759 Sin[s] 
Sin [t] +0.5606601717 Cos [t ] &&z==-2 . 828427124 Cos [s] Sin [t] +1.060660171 Sin [s] 
Sin [t] +1.224744871 Cos [t ] &&x==u+l&&y==3 u&&z==u+l , {x,y,z,s,t,u}] 

Solve: :ifun: Inverse functions are being used by Solve, so some solutions 

may not be found; use Reduce for complete solution information. More... 

{ { x->-0. 02 8 64 35 6966104 0 64 ,y->-3. 085 9307 08 983122, z->-0. 02 8 6435 6966104 064, 
u->-l. 02 8 64 35 6966104 05, t->-2 . 03 968 65050 607 42, s^l . 01 60172 80532 34 66}, 

{ x->-0. 02 8 64 35 6966104 0 64, y->-3. 085 9307 08 983122, z->-0. 02 8 64 35 6966104 0 64, 
u->-l. 02 8 64 35 69 6 61 04 05, t->2. 03 9 68 65050 607 42, s->-2. 12 557 5373057 4 4 67}, 

{x^l .28 90 64 3961420396, y^0 . 8 6719318842 61188, z->l . 2 8 90 643 96142 03 96, 
u^0 .2890643961420396, t->-0 . 121057 97 7 682 31857 , s^-1 . 72 7 302 6137 8502 34 } , 

{ x->l. 28 90 643 9614203 96, y->0. 8 67 19318842 61188, z->l. 28 90 643 9614203 96, 
u^0 .2890643961420396, t->0 . 121057 97 7 682 31857 , s->l . 4142 9003 9804 7 7 } } 

{{x,y,z}, {s,t}, {u}}/ . { sol [ [1] ] , sol [ [3] ] } 

/ (-0.0286436, -3.08593, -0.0286436} (1.01602, -2.03969} {-1.02864} \ 

l {1.28906, 0.867193, 1.28906} {-1.7273, -0.121058} {0.289064} / 

aa=ParametricPlot3D [{ 1 . 414213562 Cos[s] Sin [t] -1 . 719669914 Sin[s] 

Sin[t]+1. 478397839 Cos[t], 2.449489742 Cos[s] Sin [t] +2 . 217596759 Sin[s] 
Sin[t]+0. 5606601717 Cos [t] , -2 . 828427124 Cos[s] Sin [t] +1 . 060660171 Sin[s] 
Sin[t]+1. 224744871 Cos[t], 

FaceForm[Hue [s/Pi] ] } , { s , -Pi , Pi} , { t , -Pi , Pi} , Lighting -a False, 

Axes -a True, AxesLabel -» (StyleForm[# , "Subtitle" , FontSize^l2 , 

FontWeight->" SemiBold" ] &/0 { "X" , "Y" , "Z" } ) , 

PlotLabel->StyleForm[ "Geometry" , FontFamily->" Times" , 

FontSize^l8, FontWeight->"Bold" ] ] ; 

bb=ParametricPlot3D [ {u+1 , 3u,u+l , {Blue , AbsoluteThickness [2] } } , {u, -2 , 2 } ] ; 
cc= Show@Graphics3D [ { {AbsolutePointSize [10] , Red, Point [ {x,y , z } / . sol [ [1] ] ] } , 
{AbsolutePointSize [10] , Red, Point [ {x,y , z } / . sol [[3]]]}}]; 

Show [aa,bb, cc] 
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We have studied our lection in the first example, so we can start immediately solving the 
problem - now using immediately the Black Box (Try to do without!): 

#28: el 12 := [1 .414213562. CQS(s) *SIN(t) - 1 . 719669914. SIN(s) .SIN(t) + 1 .478397 839.CCS(t) , 

2. 44 9489742. CCS (s).SIN(t) + 2. 217596759. SIN(s) .SIN(t) + 0. 5606601717 .COS(t) , - 
2. 828427124. COSCs).SIN(t) + 1 . 060660171 .SIN(s) .SIN(t) + 1 . 224744871 -COS (t)] 


#29: line2 := [u + 1, 3>u, u + 1] 

sys2 := VECTOR(ell2 - Hne2 , k, 3) 

#30: k k 

#31: sys2 := [1 .414213562. CCS(s) *SIN(t) - 1 . 719669914. SIN(s) .SIN(t) + 1 .478397 S38.CCS(t) - u - 1, 
2. 44 9489741. CCS (s).SIN(t) + 2. 217596759. SIN(s) .SIN(t) + 0. 5606601717.C05(t) - 3*u, - 
2. 828427124. CCS(s).SIN(t) + 1 . 060660171 .SIN(s) .SIN(t) + 1 . 22474487. CCS(t) - u - 1] 


(REVERSE (NEWTONS (sys2, [u, s, t], [0, 0, 0]))) 

#32: [1, 2] 


#33: 


#34: 


' -1.028643569 -2.125575373 2.039686505 ' 

. -1.028643569 -2.125575373 2.039686505 . 

(REVERSE (NEWTONS (sys2, [u, s, t], [1, 1, 1]))) 

[ 1 , 2 ] 


#35: 

#36: 

#37: 

#33: 


' 0.2890643961 1.414290039 0.1210579776 ' 

. 0.2390643961 1.414290039 0.1210579776 . 

u_s := [-1.028643569, 0.2890643961] 
VECTOR([1-ine2], u, u_s) 

' [-0.028643569, -3.085930707, -0.028643569] 
. [1.289064396, 0.3671931332, 1.239064396] 


Mathematica Graph 


DERIVE Graph (Greyscale) 



